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Abstract 

We perform Monte Carlo simulations using the Wolff cluster algorithm of multiple q — 2, 3, 4 state 
Potts models on dynamical phi-cubed graphs of spherical topology in order to investigate the c > 1 
region of two-dimensional quantum gravity. Contrary to naive expectation we find no obvious signs 
of pathological behaviour for c > 1. We discuss the results in the light of suggestions that have been 
made for a modified DDK ansatz for c > 1. 
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1 Introduction 



In a recent paper jjj, we carried out numerical simulations of the q — 2(Ising) and q = 3,4, 10 state Potts 
models on dynamical phi-cubed graphs of spherical topology in order to see if the predictions of KPZ Q 
and DDK || for the critical exponents of the q — 2, 3,4 models were in agreement with "experiment". 
We found that the predicted and measured exponents were in quite good agreement for these models and 
that the q — 10 Potts model still displayed a first order transition on the dynamical phi-cubed graph as 
it did on a fixed lattice. We also explored the intrinsic geometry of the graphs in the simulation, finding 
a linear relation between the number of rings of length three at the critical temperature and the central 
charge, and that the peak in the ring distribution moved towards (3 C as q increased []. In the current paper 
we employ the same methods to explore multiple copies of q — 2,3,4 state Potts models on dynamical 
phi-cubed graphs. We shall simulate 2,4,8 and 16 Ising models, 5 and 10 q = 3 Potts models, and 2,4 
and 8 q — 4 Potts models all of which have, at least naively, c > 1. 

The formulae of KPZ/DDK give nonsensical results for c > 1, suggesting there is some sort of phase 
transition at c = 1. The results of simulations of dynamically triangulated random surfaces, where the 
phi-cubed graphs or their dual triangulations are embedded in d > 3 (c = d for the Polyakov action for 
a surface) dimensions seem to support this view since they give rise to very crumpled and spiky surfaces 
with no continuum limit We can include an extrinsic curvature squared or mean curvature [||, 
term in addition to the Polyakov action which may cure these ills. There have been suggestions Q that 
there is a liberation of curvature singularities at c = 1 when such a term is not included which give 
rise to the spikiness observed in simulations. As this is essentially an intrinsic phenomenon we would 
expect to see some sort of pathology for any model on a dynamical triangulation (or its dual) whose naive 
continuum limit, if it existed at all, had c > 1. The simplest sort of matter we can put on a dynamical 
graph is a spin variable, as in the various Potts models, and we can arrange c > 1 by taking multiple 
copies. The copies interact with each other only through the graph. The simplicity of the matter allows 
us to explore quite large graphs by comparison with simulations of embedded surfaces. As the graphs 
are no longer embedded the pathology, if any, cannot manifest itself as "crumpling" . The object of the 
simulations described in this paper is to see if things go wrong in some other manner for models with c 
(naively) > 1. 

The partition function we are simulating is that for multiple Potts models on dynamical random 
graphs with a fixed number of nodes N M and is given by 



z N = 53 E ex p I -|2J W I w 

where C?W is the adjacency matrix of the phi-cubed graph, S is a Kronecker delta, (3 is the inverse 
temperature 1/T and the sum over a is a sum over the different copies of the Potts model. We have, 
in general, q species of spin taking the values 0, 1, q — 1. A single Ising model has been simulated on 
both dynamical triangulations B , pj and phi-cubed graphs |llj and satisfactory agreement between the 
measured and theoretical values of the critical exponents found. The theoretical values come from both 
the work in ||,|§, as outlined below, and the exact solution If we denote the critical temperature 
for a continuous spin-ordering phase transition by T c and the reduced temperature \T — T c \/T c by t then 
the critical exponents a, (3, 7, v, 5, rj can be defined in the standard manner as t — > 

C~r a ; M~t p ,T<T c 



R N 



M(H,t = 0) ~ H 1/s , H -> 
1 

\d-2+ V ' 



<M(x)M(y)> ~ x _y, , 121 



where C is the specific heat, M is the magnetization, x is the susceptibility, £ is the correlation length and 
H is an external field. In theories without gravity it is possible to calculate a and f3 using the conformal 

1 In Jlj we stated that the peak moved towards /3 C as c increased, as c increases with q for the single Potts models. We 
shall see in this paper that this is not strictly accurate for multiple Potts models. 
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weights of the energy density operator and spin operator (for a review see |13|). Given these we can now 



use the various scaling relations 1 14 



a = 2 — vd 

= Z( d -2 + rj) 



7 = v{2-T]) 
2- 
T+?7 



to obtain the other exponents. 

Miraculously, this approach carries over directly when one considers the same theories coupled to 
gravity (i.e. on dynamical meshes): the modified conformal weights A for the energy and spin operators 
when coupled to gravity are given in (2) , Q in terms of the original conformal weights Aq by 
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A-A = -|-A(A-1), (4) 



where 

1 



£ = - — ^/25^-VT^), (5) 

and these simply replace the standard values in the calculation of the critical exponents a and (3. The 
one exact solution that is currently available, for the Ising model on a dynamical mesh, has exponents 
which satisfy the scaling relations so we assume their validity to derive the other exponents from these 
two. 

In |l| we found that there was satisfactory agreement with the results of Q , || for single copies of the 
q = 3,4 state Potts models and that the q = 10 Potts models had a first order transition. For single Potts 
models this suggests the following appealing picture: for q < 4 they display a continuous phase transition, 
which allows us to define a continuum limit, whereas the q > 4 models, as on fixed lattices, have first 
order transitions with no continuum limit and associated conformal field theory. It is also tempting to 
speculate that the difficulties that appear in the continuum formalism for c > 1 may be due to a similar 
effect, with multiple copies of the Potts models having first order transitions when c > 1. In what follows 
we shall see that the simulations fail to oblige our expectations. 



2 Spin Model Properties 

We perform a microcanonical (fixed number of nodes) Monte Carlo simulation on phi-cubed graphs 
with N = 50, 100,200,300,500, 1000 and 2000 nodes at various values of (3 between 0.1 and 1.5. The 
Monte Carlo update consists of two parts: one for the spin model and one for the graph. For the Potts 
model we use a cluster update algorithm since it suffers from much less critical slowing down than the 
standard Metropolis algorithm. There are two popular cluster algorithm implementations - Wolff jl5| 
and Swendsen-Wang |l6|. As they are equivalent for the usual two-dimensional q = 2, 3 Potts models 
we use Wolff's variant because it is computationally faster. For the graph update we use the Metropolis 
algorithm with the standard "flip" move. As we are working with phi-cubed graphs the detailed balance 
condition involves checking that the rings at either ends of the link being flipped have no links in common. 
This check also eliminates all graphs containing tadpoles or self-energies. After each Potts model update 
sweep we randomly pick NFLIP links one after another and try to flip them. After testing various values 
of NFLIP to ensure that there were enough flips to make the graph dynamical on the time scale of the 
Potts model updates we set NFLIP = N for all the simulations. 

Just as in our simulations of single Potts models we measure all the standard thermodynamic quantities 
for the spin model: energy E, specific heat C, magnetization M , susceptibility x an d correlation length 
£; and several properties of the graph: acceptance rates for flips, distribution of ring lengths and internal 
fractal dimension d. To determine v (actually vd) and (3 C separately, instead of from the usual three- 
parameter finite-size scaling fit, for example £ = £,a(\T — T c \/T c )~ v , we used Binder's cumulant jl7|j . This 
is done as follows. Binder's cumulant Un on graph with N nodes is defined as 



where < M 4 > is the average of the fourth power of the magnetization and < M 2 > is the average of 
its square. For a normal temperature-driven continuous phase transition Un - > for T > T c because 
M is gaussian distributed about at high temperature, and Un — > § for T < T c because a spontaneous 
magnetization M sp develops in the low temperature phase. At T — T c , U n has a non-trivial value which 
scales with N according to 

Un — tN^<* . (7) 

Therefore the slope of Un with respect to T (or f3) at T c gives This is not much use as it stands 
since it involves knowledge of T c , but the maximum value of the slope scales in the same way, so we can 
extract vd from 

, dU /v , i , * 

max(^)^iV-. (8) 

We show the results of doing this in Figs. la,b,c for the q — 2,3,4 state Potts models respectively. From 
the figures it is clear that there is little variation in the values of vd for a given q as the number of copies 
of the spin model is increased. We therefore use a weighted average vd for each q in what follows: 3.6(4) 
for q = 2, 2.9(4) for q = 3, and 2.8(2) for q = 4. 

We now make use of these values of vd in the standard finite-size scaling relation (with L replaced by 
Nd since we do not know d a priori) 

|/^-/?f|~JV-^ (9) 

to extract , using obtained from the position of the maximum in the slope of Un or from the 
peak in the specific heat C. The end results are shown in Table 1, where we list the model in the form 
'number of copies' q-value, central charge c, vd, (3 C estimated from Un and f3 c estimated from C. It is 
clear that the peak in the specific heat gives the more accurate estimates for (3 C , and we use this in the 
finite size scaling analysis which follows. 



model 


c 


vd 


P?{Un) 


PT(C) 


2 q2 


1 


3.5(3) 


0.80(5) 


0.79(1) 


4 q2 


2 


3.9(4) 


0.87(4) 


0.82(1) 


8 q2 


4 


3.8(4) 


0.84(4) 


0.84(1) 


16 q2 


8 


3.4(5) 


0.92(5) 


0.88(1) 


5 q3 


4 


2.9(2) 


0.93(3) 


0.91(1) 


10 q3 


8 


3.2(5) 


0.99(3) 


0.96(1) 


2 q4 


2 


2.8(2) 


0.96(3) 


0.95(1) 


4 q4 


4 


2.9(2) 


0.97(3) 


0.96(1) 


8 q4 


8 


2.8(1) 


1.02(3) 


0.99(1) 



Table 1: Fitted values of vd and inverse critical temperature /3 C 
from Binder's cumulant Un and specific heat C . 

We have chosen this slightly non-standard approach to obtaining f3 c , rather than looking for the 
crossing of cumulants on different sized meshes, because we found it very difficult to judge the crossing 
points accurately. Applying our approach to the single Ising model in [|lj gave an accurate estimate of 
its known j3 c . A comparison of the results of multiple Potts models with those for single models in JjJ 
shows that C increases slowly with the number of models, but that vd stays roughly constant at a value 
somewhat larger than that for a single model. 

We can now measure the other critical exponents either from the singular behaviour of the thermo- 
dynamic functions 

C = B + C t- a , M = Mot?, X = Xot-\ £ = (10) 
knowing (3 C , or from the finite-size scaling relations 

C = B' + C' Q N^, M = MqN^t, x = x' N^ (11) 

using the previously obtained values of vd. The results obtained for [3/vd, ^ /vd, 7 and v in the various 
models obtained from the above fits are shown in Table 2 
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model 


c 


(3/vd 


"f/vd 


7 


V 


d 


2 q2 


1 


0.11(1) 


0.83(1) 


2.4(1) 


0.87(2) 


2.33(4) 


4 q2 


2 


0.08(1) 


0.94(2) 


3.1(2) 


0.93(4) 


2.75(4) 


8 q2 


4 


0.09(1) 


0.98(1) 


3.3(2) 


0.83(2) 


2.67(4) 


16 q2 


8 


0.08(1) 


1.05(2) 


2.7(1) 


1.02(2) 


2.59(4) 


5 q3 


4 


0.09(1) 


0.91(1) 


2.2(1) 


0.97(2) 


2.58(3) 


10 qo 


o 
8 


0.05(1 ) 


1.01(2) 


2.8(2 ) 


1.24(2) 


2.56(3) 


2 q4 


2 


0.09(1) 


0.86(1) 


2.0(1) 


1.16(3) 


2.63(3) 


4 q4 


4 


0.11(1) 


0.90(1) 


1.7(1) 


0.85(2) 


2.52(3) 


8 q4 


8 


0.16(1) 


0.98(1) 


2.9(1) 


0.85(4) 


2.47(3) 



Table 2: Measured values of critical exponents (3^,v 
and internal fractal dimension d from N — 2000 graphs. 

We see immediately that 'y/vd increases for q = 2,3,4 with the number of models, whereas (3/vd 
decreases for q = 2, 3 and increases for q = 4. However, the measurements of (3/vd proved to vary the 
most from the predicted values of KPZ for the single potts models in [Q, so one should really conduct 
simulations on larger meshes before taking the trends in [3/vd at face value. In order to make a meaningful 
comparison of these results with those for single Potts models in pj, where j/vd decreases with q and 
(3/vd increases with q, we should look at multiple models with the same c. Then we see that for c = 2 
j/vd decreases from 0.94(2) for q = 2 to 0.86(1) for q = 4, and (3/vd increases from 0.08(1) for q = 2 
to 0.09(1) for q = 4. Similarly, for c = 4 as q increases from 2 to 3 to 4 j/vd decreases from 0.98(1) 
to 0.91(1) to 0.90(1), and (3/vd increases from 0.09(1) to 0.09(1) to 0.11(1). And finally for c = 8 -y/vd 
decreases from 1.05(2) to 1.01(2) to 0.98(1), and (3/vd increases from 0.08(1) to 0.05(1) to 0.16(1). From 
these trends we would argue that the value of (3/vd for 10 q3 is a little low, but otherwise things look 
consistent. 

The direct fits to 7 are poorer than those to (3/vd and 7 /vd due to the uncertainty in (3 C . In fact, 
if we increase (3 C from the values given in the last column of Table 1 by the size of the error bar (0.1) 
and re-do the fits we obtain values for 7 which are 0.2 lower than those given in Table 2. Therefore our 
values for 7 can be trusted only to about 10%. These values are all larger than those obtained for single 
Potts models - which makes sense as 7 controls the divergence of the susceptibility and obviously with 
more spin models this divergence will be stronger. For single Potts models |lj we know that 7 decreases 
with q so as above we should compare our multiple models with the same c. For c = 2 and c = 4 we find 
decreasing values but for c = 8 we have 7 = 2.7,2.8,2.9 for q ~ 2,3,4 respectively. We could say that 
this casts doubt on the fits for 10 q3 and 8 q4 but given the errors judgement should be reserved until 
the next generation of simulations is performed. As a consistency check we can multiply our fit values of 
y/vd by vd and see how close we get to the fitted values for 7. Surprisingly, almost all results are within 
10%, and the two worst (16 q2 and 4 q4) only differ by 30%. 

We have not attempted to fit the exponent (3 directly here since even for single Potts models we never 
managed to obtain satisfactory results for this [Q. 

We fitted the exponent v from the power law divergence of the correlation length £ at (3 C . However, 
£ itself must be obtained from a fit: the 2-point correlation function T should behave as 

r(r) = 5> if 7 i+r = ce~ mr , (12) 

where m = l/£, and the sum is over some number of measurements made on each graph with the position 
of the spin Ui being chosen randomly, r is the internal distance between two spins on the graph, i.e. the 
fewest links between them. As two fits are involved, both requiring knowledge of (3 C , the results are not 
particularly reliable. If we re-do the fits with (3 C changed to (3 C + 0.1 then the sixth column of Table 2 
would read instead 1.01(2), 1.11(5), 1.00(2), 1.19(2); 1.15(3), 1.55(2); 1.45(4), 1.04(2), 1.10(2). Combining 
both sets of fits (we do not consider fits using (3 C — 0.01 since our other estimates of (3 C from Binder's 
cumulant are always larger than those from the specific heat) yields our best estimates for these v values 
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as 1.0(1), 1.0(1), 0.9(1), 1.1(1); 1.0(1), 1.4(1); 1.3(1), 1-0(1), 1.0(1). Thus we conclude that v is 1 to an 
accuracy of 10% for most if not all the models. 

As a further consistency check one can also measure the internal fractal dimension d directly, as we 
discuss in the next section, and see how the product of the (less accurate) individual measurements of v 
and d compares with the direct determination of vd from Binder's cumulant. We can see by comparing 
Tables 1 and 2 that the product underestimates vd for almost all the models, but still appears to be fairly 
good considering the multiple fits involved. 

Unfortunately, due to the extra adjustable constants (B, B') in eqs. [H],[ll] for the specific heat, we 
are unable to make unambiguous fits for a and ajvd. Although we cannot fit directly for a we can use 
our values of vd and assume that the scaling relation a = 2 — vd is still valid, which gives (with the 
averaged value of vd for a given q) a — —1.6(4) for q = 2, a — —0.9(4) for q = 3 and a = —0.8(2) 
for q = 4. We thus see that this gives third order transitions for all q, which are becoming weaker with 
increasing q. These results should be compared with the single Potts models where we have (analytically) 
a = — 1 for q = 2, a = —0.5 for q — 3 and a = for q = 4. Adding multiple copies at a given q thus 
appears to make the third order transitions stronger for q = 2, 3 and change the borderline q = 4 case 
from logarithmic second order to third order. In addition, visual inspection of the graphs of specific 
heat for various N suggests that they are not diverging with N, as one might expect for a higher than 
second order transition. We can, of course, rule out a first order transition because we do not observe 
any discontinuities in the energy or magnetization. 



3 Graph Properties 

We are particularly interested in the properties of the graphs when c > 1 to see if there are any indications 
of pathological behaviour as the spin models themselves, by the evidence of the preceding section, seem 
remarkably well behaved. We measure the acceptance rate for the Metropolis flip move to confirm that 
our graphs are really dynamical. The flip can be forbidden cither from the graph constraints coming 
from the detailed balance condition or from the energy change of the spin model, so we can decompose 
the flip acceptance rate into two parts: AL - the fraction of randomly selected links which can be flipped 
satisfying the graph constraints; and AF - the fraction of links satisfying the graph constraints which 
are actually flipped, i.e. pass the Metropolis test using the Potts model energy change. AF and AL 
are shown for the 2000-node simulations versus the reduced temperature t in Figs. 2 and 3 respectively. 
These are interesting because they show that the dynamical properties of the phi-cubed graphs are 
apparently determined by the central charge of the models that live on them, even for values of the 
reduced temperature t some distance from the critical point. The various spin models group into three 
separate curves determined by the value of their central charges, c = 2, 4, 8. The characteristic dip in the 
flip acceptances away from the phase transition that is present in the single Potts models and dynamically 
triangulated random surface models near the crumpling transition is clearly visible. For a given q the 
position of the dips stay roughly fixed as the number of models is increased, though they do appear to 
edge slightly closer to t = as q is increased. 

We can also examine the distribution of ring lengths in the graph, which is the discrete equivalent 
of measuring the distribution of local Gaussian curvatures in the continuum. For pure quantum gravity 
(no spin model living on the graph) it is possible to analytically calculate this jl8) . The probability P of 
finding a ring of length I is given by 

P^ oo (0 = 16(-) lKl _ 1){ (13) 

which decays exponentially as I increases. The minimum possible ring length is 3. If we plot the fraction 
of rings of length three (PR3) in Fig. 4 and compare with AF and AL (Figs. 2 and 3) we see that it has 
a peak very close to the dip in AL. This is reasonable since both PR3 and AL depend only on the graph, 
whereas AF depends on the Potts model. The suggestion made in Q, namely that the peak in PR3 
approached [3 C from below as c — * 1, is not therefore bourne out by our further simulations of multiple 
models. Just as for the flip acceptances, the (3 values of the peaks stay roughly fixed for a given q as 
the number of models is increased, though they get slightly closer to t — with increasing q. It is also 
obvious from Fig. 4 that the PR3 curves split into three again for c = 2, 4, 8 for all the plotted values of 
t, in spite of the small shift in the peak positions with q. 
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It is also interesting to examine the height of the peaks in PR3 (i.e. maximum values of PR3) and the 
values of PR3 at the critical points as a function of c. If we plot the difference in the fraction of rings of 
length three at the critical point (3 C from the pure gravity fraction (eq. [l^ with I = 3) against the central 
charge we find a straight line with slope 0.004(1) and intercept 0.004(1). We can also plot a difference 
using the peak height in the fraction of rings of length three against the central charge to obtain another 
straight line with slope 0.014(1) which passes through the origin. This peak does not occur at j3 c so 
the correlation length is not infinite and we cannot expect the results of conformal field theory to apply 
there, but a linear relation with c persists. The fits are shown with the data in Fig. 5. This should be 
compared with a plot of similar data for smaller c values in Fig. 7 of Q, which agrees with the data here 
in that range. 

We complete our discussion of the graph properties by measuring their internal fractal dimension d. 
We use the most naive definition of distance (the fewest links between two nodes) so we are considering the 
"mathematical geometry" rather than the "physical geometry" in the terminology of |Q 0. We measured 
the internal fractal dimension at the critical values of j3 for all of our models on the N = 500, 1000, 2000 
node graphs and then extrapolated the results to graphs of infinite size obtaining the values listed in 
the last column of Table 2. From these we can see that the internal fractal dimension varies little, if at 
all, with c and q and all are fairly close to the value of approximately 2.8 that we found in [Q for single 
Potts models. It thus appears that although the curvature distribution may be becoming more singular 
(slowly) with increasing c, as is witnessed by the increase in PR3, the fractal dimension is changing very 
little. 

4 Conclusions 

We have simulated multiple copies of q = 2, 3, 4 state Potts models on dynamical phi-cubed graphs all of 
which have, at least naively, c > 1. We have no evidence for a sudden change in behaviour in the spin 
models as c is increased through 1, but instead found only a gradual change in the critical exponents 
we have measured. The idea that a first order transition could explain the breakdown of the KPZ/DDK 
results for multiple Potts models with c > 1, as it does for q > 4 Potts models does not appear to be 
supported by our simulations in which the multiple Potts copies behave much as their single brethren do. 
That the spin models do interact with each other is clear from the change in the exponents and critical 
temperatures from those of the single models. The properties of the phi-cubed graphs such as the flip 
acceptances and the fraction of rings of length three which vary with the number of models also bear this 
out, and indicate that the central charge alone determines the internal geometry, even for c > 1. 

It is possible that the failure of the KPZ/DDK results for c > 1 is merely a breakdown of the ansatz 
rather than an indication of a completely new phase. Suggestions to this effect have been made in pl[ 
where an extended ansatz involving the coupling of the diffcomorphism ghosts to the Liouville theory was 
suggested. It would be most interesting to explore some of the suggestions made for pinning down the 
parameters of the ansatz in j2l[] to see if it is possible to calculate any of the exponents measured here 
and to reproduce their slight variation with c. On the numerical side it would obviously be interesting to 
extend the work in the current paper to larger meshes and better statistics, in particular to attempt to 
obtain values for a or ajvd directly, rather than using the scaling relation. It would also be a worthwhile 
exercise to measure the string susceptibility 7 s t r for our multiple Potts models, as this is one of the easiest 
exponents to calculate analytically. 
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Figure Captions 

Fig. 1. Fits to maximum slope of derivative of Binder's cumulant versus N to extract vd for a) q=2, b) 
q=3 and c) q=4. 

Fig. 2. Flip acceptance AF as function of reduced temperature t for all q. 
Fig. 3. Flip acceptance AL as function of reduced temperature t for all q. 

Fig. 4. Probabilities of rings of length three PR3 as function of reduced temperature t for all q. 

Fig. 5. Difference in PR3 at /3 C , and at its peak, from the pure quantum gravity value versus the central 
charge c, for multiple Potts models. 



8 



